Soft Modes at the Stacking Faults in SiC Crystals 
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We use first-principles calculations based on density functional theory to determine and under- 
stand the driving force of the observed stacking fault expansion in SiC. We verify the suggestion 
based on recent experiments that the free energy difference between the faulted and the perfect 
structures is responsible for this expansion, and show that its origin lies in a large entropy associ- 
ated with soft vibrational modes of the faulted SiC structure. These soft modes involve shearing of 
SiC on a long length-scale and are absent in related group IV semiconductors, such as Si, Ge and 
C. 



Among the group IV semiconductors (Si, Ge, C, SiC), 
SiC is unique in exhibiting hundreds of polytypes that 
differ only in the stacking sequence of Si-C double layers^. 
Due to its excellent physical properties such as high hard- 
ness, low density and thermal expansion, it has found 
numerous applications in automotive components, seal 
faces, armour, mechanical seals, bearings, heat exchanger 
tubes, etc^^ As it is available in large single crystal wafer 
form^, SiC is also interesting as a wide band-gap semi- 
conductor for electronic applications: it is a promising 
material for the next generation power diodes owing to 
its low on-resistance, high break-down voltage, radiation 
stability and low switching losses^. Recently, there have 
been efforts aimed at developing dilute magnetic semi- 
conductors (DMS)^ based on SiC for use in spintronic 
devices. 

However, formation of line defects^^^^ and expansion of 
stacking faults^ in SiC crystals are known to be possible 
causes of degradation of the performance of SiC-based 
power diodes. Understanding the mechanism of stacking 
fault expansion in SiC is important to its use in electronic 
devices and it may also shed light on a related fundamen- 
tal question of the stability of polytypes of SiC^. It is also 
interesting to find out if this mechanism operates in other 
group IV semiconductors. 

In a power electronic device current densities as high 
as 50-100 A/ cm? lead to an enormous amount of heating 
and operational temperatures are higher than the normal 
ones. Optical emission microscopy (OEM) based experi- 
ments on 4H-SiC used in high current density conditions 
revealed stacking fault expansion in the basal plane of the 
system, establishing that hexagonal polytypes of SiC are 
unstable at high temperatures with respect to formation 
and expansion of stacking faults^. Based on this work, 
free energy differences was suggested to be the driving 
force of stacking fault expansion. 

Many groups have carried out first-principles investiga- 
tion of the stability of various polytypes of SiC and their 
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stacking fault^iiiii^. With very small energy differences 
involved, the exact ordering energies of polytypes is a bit 
scattered in these studies. However, they all estimated 
very small but positive stacking fault energies. It is not 
clear yet what drives the stacking fault expansion in SiC 
at finite temperature. 

In this work, we use first-principles calculations to de- 
termine free energy differences between different poly- 
types and stacking-faulted structures of SiC through de- 
termination of configurational and vibrational entropies. 
Configurational and vibrational contributions to free en- 
ergy are determined from the knowledge of full stacking 
fault energy surface^^ and of the vibrational spectrum re- 
spectively. We show that soft shear modes developed in 
the faulted structures drive their stability at finite tem- 
perature. 

Diamond (or zincblend, labeled here as 3C) structure, 
the common polytype of all group IV semiconductors, can 
be described as two interpenetrating FCC lattices sepa- 
rated by [lll]a/4, a being the lattice constant. Along 
[111] direction, triangular lattices of atoms described 
with high symmetry positions A, B and C are stacked pe- 
riodically as (AaBbCc), where upper and lower case let- 
ters label the two FCC sublattices. In case of SiC, carbon 
atoms occupy sites of the second sublattice. Atoms in the 
polytypes obtained with a variation in this stacking se- 
quence (studied in this work) are all tetrahedrally coordi- 
nated. In the Hagg notation^^ used commonly in the rep- 
resentation of stacking sequences, (AaBb), (BbCc) and 
(CcAa) are denoted as (+), and the anti-cyclic stacking 
sequences eg. (BbAa) are denoted as ( — ). Thus, 3C- 
SiC (with AaBbCc stacking), 4H-SiC (with AaBbAaCc 
stacking) and 6H-SiC (with AaBbCcAaCcBb stacking) 
polytypes are represented with periodic units of (+ + +), 

(+ H ) and (+ + H ), respectively. In the 

Hagg notation, a stacking fault is introduced simply 
through a change in sign. Deformation type stacking 
faults in 3C-SiC, 4H-SiC and 6H-SiC systems, which are 
bordered by Shockley partials on the basal plane, are 

(...+++!-+++++...), (...++--!- + --++--...) 

and (...+ + H + H ^ + H --O^ respec- 
tively, where the symbol '|' indicates the stacking fault 
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plane(lO). There are two types of (111)/(0001) planes 
where a deformation type stacking fault can form: (a) 
A\a called as "shuffle", and (b) a\B called as "glide". 

Our calculations are based on first-principles pseu- 
dopotentials within Density Functional Theory as im- 
plemented in the Plane Wave Self Consistent Field 
(PWSCF)^^ code. We employ a local density approxi- 
mation (LDA) to exchange correlation energy functional 
and use ultra-soft pseudopotentials^^ to represent inter- 
action between ionic cores and valence electrons, and a 
plane wave basis with an energy cut off of 30 Ry (180 
Ry) in the representation of the Kohn-Sham wavefunc- 
tions (density). Supercells with 12 atomic planes for ele- 
mental semiconductors and 24 atomic planes for SiC are 
used in the calculations of configurations with of stacking 
faults. Corresponding Brillouin zones are rather small in 
the direction (perpendicular to the plane of the fault) 
and integrals over them were sampled with a uniform 
5x5x1 mesh. Positions of atoms inside the supercell 
were relaxed to attain a minimum energy structure us- 
ing Hellman-Feynman forces in the Broyden, Fletcher, 
Goldfarb, Shanno (BFGS)-based method^^. Structural 
parameters for Si, Ge and C in diamond structure and 
3C, 4H and 6H polytypes of SiC agree with the experi- 
mental values within the typical LDA errors. We find the 
cohesive energies of all the SiC polytypes are rather close 
to each other and both 6H and 4H structures are almost 
equal within the calculational errors and the lowest in 
energy. 

We first determined generalized stacking fault energy 
(7—) surfaces^ for the (111) or basal planes of all the 
crystals studied here. It corresponds to the energy re- 
quired to displace a half of the crystal on one side of the 
given plane relative to the other and exhibits periodicity 
of the crystal plane. For a given displacement path go- 
ing from the origin to a burgers vector b, a maximum in 
the 7— surface gives the the energy barrier for slip along 
that path. The lowest of these barriers among the family 
of various paths joining origin and b is called the un- 
stable stacking fault energy^^ 7^^^. A local minimum of 
the 7— surface defines an intrinsic ( also called deforma- 
tion(lO) ) stacking fault. We used fourier transform to 
analyze the 7— surface by sampling it on a uniform 5x5 
mesh in the planar unit cell. 

In all the tetrahedral systems considered here, our re- 
sults rule out the possible slip on shuffle plane as it is 
energetically 25-36 times more expensive than the slip 
on the glide plane. The 7— surfaces of the glide plane 
of all tetrahedral semiconductors are topologically simi- 
lar (see Fig. 1 for the 7— surface of SiC). We find 7^^^ is 
about the same for the three polytypes of SiC, which is 
about half the 7^^^ of diamond and larger by 80 % than 
that of Si and Ge (see Table I). The trend in our esti- 
mates of Jus suggests that nucleation of dislocations on 
(111) planes is much easier in Si and Ge than in SiC and 
diamond^^. The energy of intrinsic stacking fault located 
at (2/3,1/3) in the 7— surface ijisf)-, bears a rather dif- 
ferent trend: it is smallest for the polytypes of SiC and 



TABLE I: Calculated and experimental estimates of intrinsic 
(jisf) and unstable (jus) stacking fault energies for the slip 
in glide plane. 



System 


7is/(calc.) 


7is/(expt.) 


Jus 




(mJ/m2) 


(mJ/m^) 


(J/m^) 


Si 


46.9 


69^^ 


1.7 


C 


250 




5.5 


Ge 


48.5 




1.6 


3C-SiC 


10.1 




2.8 


4H-SiC 


9.1 


14.7±2.5ii 


2.9 


6H-SiC 


2.6 


2.5it0.9^2 


2.9 



5 times larger for Si and even larger for diamond. These 
values generally agree well with experimental estimates 
and other calculations^^ wherever available. Significantly 
lower values of 7^5/ (and somewhat large jus) of SiC 
support large areas of stacking faults which sometimes 
extend right upto the crystal boundaries with partials 
disappearing at the surface. 

We estimated the free energies of perfect and faulted 
SiC with F = E -\- Fconfig + Fvih^ where E is the total 
energy obtained with DFT calculations, Fconfig and FyH) 
are the configurational and vibrational contributions to 
free energy respectively. Fconfig is obtained from the de- 
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FIG. 1: Contour plot of Generalized Stacking Fault Energy (7) 
Surface for slip on the glide plane of 4H-SiC 

tailed structure of basins of the 7— surface centered at 
(0,0) and (2/3,1/3) corresponding to perfect and faulted 
structures. We approximated the 7— surface in the neigh- 
bourhood of the centres (xq, ^0) of these basins (sampled 
by the system at finite temperature T) with a symmetry 
invariant parabolic ioTia:E{x^y) = Eo -\- K{{x — XoY + 
{y — Vo)'^)^ to determine configurational contribution to 
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free energy^ 



(a) 



conf 



-keTLogi 



K 



2nkBT 



(1) 



where x and y are the components of displacements for 
which the GSF has been obtained, the parameter K was 
determined from DFT calculations carried out on a fine 
uniform 3x3 mesh over an area of 0.174 bohr^ centered 
in each basin. Similarly, vibrational free energy is deter- 
mined with 



Fvib = -ksT Log{2sinh{ 



Tiuoj 



zq 



2hT 



)) 



(2) 



coiq being the frequency of a phonon i with wave vector 
(7, obtained with DFT linear response calculations^^. We 
sampled g at F and K points with weights of | and | 
respectively, omitting acoustic modes at the F. 

The stacking fault energy as a function of temperature 
is obtained using: 7,(T) = 7,(0) + AFcon/(T) + AF^i5(T) 
where 7s (0) is the stacking fault energy at zero temper- 
ature, and AF is the difference in free energy of the 
perfect and the faulted structures. In the approxima- 
tion that phonon frequencies do not change with tem- 
perature, free energies vary linearly with temperature at 
low T For 4H-SiC, we find the vibrational contribution 
to AF/T (-0.2037 x 10"^ J/m^iT) dominating over the 
configurational contribution {0.27xl0~^ J/ m^K). Evi- 
dently, vibrational entropy plays a crucial role in stabi- 
lizing the faulted 4H-SiC structure at high temperatures. 
The stacking fault energy as a function of temperature 
(see Fig. 3) reveals that the faulted structure is energeti- 
cally favorable at 60 K, supporting the observed stacking 
fault expansion and the suggestion that it is driven by 
the free energy difference^. 

We determine the precise microscopic mechanism that 
stabilizes the faulted structure by examining mode-by- 
mode contribution to the vibrational free energies. In 
comparison with the perfect 4H structure, there emerge 
four soft phonon modes (see Figure 2) at F point in 
the faulted 411 structure with frequencies 107.3, 107.4, 
109.9, 110.0 cm~^, which lowers its free energy through 
its contribution to the entropy. We note that most of 
these soft modes involve shear straining of the structure 
over a longer length-scale in a direction perpendicular 
to the plane of the fault. These soft modes couple with 
shear strain and can drastically reduce shear elastic mod- 
uli which in this case is C44. Thus, large shear strains 
can be induced even with small shear stresses in the sys- 
tem. We note that our estimate of the temperature above 
which the faults get stabilized can have large errors due 
to many approximations made in our analysis. However, 
the mechanism we have identified here is expected to be 
reliable. In contrast to the 4H structure, we find the con- 
figurational contribution to free energy alone is adequate 
in stabilizing faults in the 611 structure. This means the 
energy basin of the perfect 611 structure is stiffer than 
the energy basin of its faulted form yielding larger con- 
figurational entropy to the faulted form. 
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FIG. 2: (a) Phonon Density of States of perfect and faulted 
4H-SiC (b) Soft shear phonon mode of frequency 107.3 cm~^ 
in faulted 4H-SiC 



Finally, we would like comment on the relative stabil- 
ity of 411 and 6H structures of SiC. Many gionp^^^^ 
have reported energies of these structures within a couple 
of meV per atom of each other. However, these are very 
small energies and we think both structures are equally 
stable at T=0 K within our calculational errors. By in- 
cluding vibrational contribution to free energy at nonzero 
temperature, we find that the 611 structure is more sta- 
ble than the 411 structure by a sizeable difference in free 
energy (with a rate of 6xl0~^ eV/atom/K). Along with 
our earlier conclusion that the 6H structure would also 
tend to form faults at finite temperature, we conclude 
that both 4H and 6H structures of SiC are metastable at 
finite temperature. 

In summary, we have verified the speculation based on 
experimental work^ that free energy difference is the driv- 
ing force of stacking fault expansion in SiC. The stability 
of the faults has been shown to arise from the vibrational 
entropy and particularly soft phonon modes that involve 
atomic displacements that shear the faulted structure on 
long length-scale. As a result, stacking faults are ex- 
pected to grow at finite temperature in both 4H and 6H 
polytypes of SiC and this is a fundamental limitation of 
SiC for use in devices. Our work should stimulate further 
theoretical and experimental work that focus on destabi- 
lizing stacking faults in SiC through suitable doping or 
other means. 
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FIG. 3: Stacking Fault energy (jisf) of 4H-SiC as a function 
of temperature 
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